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Abstract. We analyze the stress induced and thermally assisted roughening of a forest of 
superclimbing dislocations in a Peierls potential in the presence of 3 He impurities and randomly 
frozen in static stresses. It is shown that the temperature of the dip Td in the flow rate observed 
by Ray and Hallock (Phys.Rev. Lett. 105, 145301 (2010)) is determined by the energy of 3 He 
activation from dislocation core. However, it is suppressed by, essentially, the logarithm of the 
impurity fraction. The width of the dip is determined by inhomogeneous fluctuations of the 
stresses and is shown to be much smaller than Td. 



1. Introduction 

Strong interest in the supersolid state of matter in free space pQ has been revived by the 
recent discovery of the torsional oscillator (TO) anomaly in solid He [SJ. While finding no 
supersolidity in the ideal 4 He crystal, ab initio quantum Monte Carlo simulations did find 
that some grain boundaries [3], dislocations [H [5] or crystal boundaries [6] support low-d 
superfluidity spatially modulated by the surrounding lattice. In principle, a percolating network 
of superfluid dislocations [7J could explain the TO anomaly if the dislocation density is 3-4 
orders of magnitude higher than it is expected to be in a slowly grown and well annealed 
crystal. Consistent with such expectation is also a very small rate of the critical superflow 
through solid 4 He (occuring presumably along dislocations with superfluid cores) observed in 
the UMASS-Sandwich experiments [H [9]. Thus the nature of the TO anomaly in solid 4 He 
remains unclear. 

In the present work we focus on the very unexpected feature of the UMASS-Sandwich 
experiment [9] - the strong suppression of the supercritical flow rate V cr and then its recovery 
in a narrow range of temperatures. Such a feature occurs well below (about 10 times) the flow 
onset temperature To ~ 0.5 — 0.6K |9j. In our recent work [10] it has been proposed that this 
feature can be due to the stress induced roughening of superclimbing dislocations supporting the 
superflow. The key question, though, is why the temperature for the anomaly is pretty much 
independent of dislocation density and why the dip is so narrow despite disordered nature of the 
solid 4 He. Here we are proposing the explanation within the scenario where 3 He impurities as 
well as internal stresses biasing superclimbing dislocations turn out to be crucial. We will show 
that the dip-anomaly temperature is determined by the 3 He activation energy E a ~ 0.5 — 0.8K 
which is reduced by a large logarithm of 3 He fraction. 



2. 3 He precipitation on dislocation core 

Atoms of 3 He provide pinning centers for dislocations. Thus, Frank's forest of super climbing 
dislocations can be viewed as consisting of free segments of quantum superclimbing strings [5|I10] 
of length L equal to the average distance between 3 He impurities absorbed on dislocations, 
provided L is less than a typical distance between the network cross- links. Here we derive the 
equation for L as a function of temperature within the thermodynamical model including the 
total number of 3 He atoms in the crystal bulk N b and the total number of the impurities N d 
absorbed on dislocations. It is natural to presume that cores of the dislocations constitute a 
small fraction x d « 1 of the total sample. Free energy of such system can be written in the 
following form 

G = E a N b + TN b ln(N b /eN ob ) + TN d ln{N d /eN od ), N b + N d = N 3 , (1) 

where N3 is the total number of 3 He atoms (with their fraction X3 = Ns/NQ b << 1); Nqj, 
denotes the total number of sites available for 3 He atoms in the bulk; No d = x d No b stands for 
the total number of sites available for 3 He atoms on dislocations. In this model we presume 
that No d > N3 so that there is no need to consider the effect of the entropy reduction on 
dislocations when all 3 He condense on dislocations. Including this effect, while introducing 
additional technical difficulties, does not alter the main conclusion of this work. 

Varying G with respect to N d as dG/dN d = 0, while keeping N3 fixed, and introducing the 
fractional 3 He concentration on dislocations X = L^ 1 = N d /No d , we find 



x d + exv(-E a /Ty 

where the inter-impurity distance L along a dislocation core is measured in units of a typical 
inter-atomic distance b ~ 3 — 4A in solid 4 He. In what follows we will be using this value 
L = L(T) as a length of dislocation segments undergoing the stress induced roughening [TO] . 



3. Stress induced roughening of a dislocation 

Superclimbing dislocation is modeled as a quantum string oriented along the rr-axis and strongly 
pinned at its both ends x = 0, L [TTJ. Here we will be considering L determined by the mean 
distance between 3 He impurities absorbed on dislocations, Eq.([2]). The string displacement 
y(x, t) along the y-axis depends on the time t and is measured in units of the inter-atomic 
spacing (w Burger's vector b) with respect to its equilibrium y = (no tilting is considered). 
The Peierls potential induced by the crystal is taken as Up = —up cos (2-Ky(x, t)). The partition 
function Z has the form [5l \12\ [TO] 

Z = J Dy(x,t) Dp(x,t)D<f)(x,t)exp(-S), (3) 

S = / dt]T[i(p + no)V^ + ^(V^) 2 

+ ^P~yf + ^ (V t y? + V>(V x y?) 

- up cos (27ty(x,t)) -Fy(x,t)}, (4) 



where all the variables are periodic in the imaginary time t > with the period /3 = l/T 
(units h = 1, Kp = 1); the core density p and the superfluid phase (ft are canonically conjugate 
variables, with p' = p — y being the local superfluid density; the derivatives Vt, x y, ^t,x4> are 
understood as finite differences in the discretized space-time lattice (with 200 time slices and 



x = 0, 1, 2, L in units of 6), with Vt,x4> defined modulo 2ir (in order to take into account phase- 
slips); no,/?o stand for the average filling factor ( we choose no = 1) and the bare superfluid 
stiffness, respectively, with the bare speed of first sound taken as unity. 

The first two terms in Eq.Q describe the superfluid response of the core [5J, and the third 
term accounts for the superclimb effect jSJ - building the dislocation edge so that the core climb 
y — >■ y ± 1, ±2, ... becomes possible by delivering matter p — > p± 1, ±2, respectively, along the 
core [5]. The dislocation is assumed to be attached to large superfluid reservoirs at both ends, 
with spatially periodic boundary conditions for the supercurrent. 

The terms oc m in Eq.@ account for the elastic response of the string, with m and standing 
for the effective mass of the dislocation core (per b) and the bare speed of sound, respectively. 
Since the main source of kinetic energy are supercurrents, we have left out the term ~ (Vty) 2 in 
Eq. ([!]). The parameter m in Eq.(j3|) is not actually a constant. It contains a contribution from 
the Coulomb-type interaction potential oc l/\x\ between jogs (or kinks, cf. [13]) separated by a 
distance x [13] • Accordingly, m has a logarithmic divergent factor with respect to a wave vector 
q along the core m(q) = ■ 1 + Uc In (\ + jj^yi^j , where mo is of the order of the atomic 

4 He mass and Uc ~ 1 is a parameter characterizing the strength of the interaction [131 [12] . In 
solid 4 He, the zero-point fluctuation parameter K = irh / (AniQbVd) ~ 1 [13] . which justifies the 
necessity of implementing full quantum-thermal analysis. We present our numerical results for 
Uc = 1>^ 2 = 5, K = 1. It is important to note that the main results are not qualitatively 
sensitive to the long-range interaction. 

As discussed in Ref.|10]. in the case of a single dislocation in an ideal crystal the linear force 
density F ~ ba (ignoring spatial indices) is determined by applied chemical potential in the 
setup [9j. Here we consider F in Eq.(|3]) determined by external stresses a ~ 0.1 — O.Olbar 
usually existing in solid 4 He ( a = F in units of b). Such stresses are inhomogeneous. We argue 
that they are responsible for the roughening effect [10] in real samples, and their inhomogeneity 
may actually wash out the effect of the applied chemical potential. 

Monte Carlo simulations have been conducted with the Worm Algorithm (WA) [15] for the 
superfluid part of the action, with the Peierls term treated within the Villain approximation 
similarly to Refs. \13\ I12j. The renormalized superfluid stiffness p s (T,F) and compressibility 
k(T,F) have been calculated in terms of the windings of the dual variables |15[ [T3"l 112] . No 
significant effect of the bias F was found on p s (T, F), and thus we consider p s (T, F) = p s (T, 0) = 
Ps(T). In contrast, k(T,F) does experience quite dramatic renormalization as found in Ref.|10] 
and demonstrated in Fig. [TjThus,the renormalized speed of first sound 



exhibits the dip |10] , According to the Landau criterion such speed limits the critical superflow 
rate along the dislocation core. 

The bare stiffness po{T) vanishes above some temperature To comparable to the bulk A- 
temperature (~ 1 — 2K). We have used To = 0.2 (in the dimensionless units as a fraction of 
the Debye temperature To for the first sound), and considered low temperatures - such that 
p s (T) stayed unchanged within 1-10% of its T = value. In other words, the thermal length 
Lt ~ Ps{0)/T, above which p s {T) becomes suppressed, is the largest scale in the problem. 
Jog-antijog pairs as quantum objects can be created spontaneously by a macro scopically small 
stress F > F c oc 1/L applied to a superclimbing dislocation of length L - analogous to the 
creation of kink-antikink pairs along a stressed gliding dislocation [14} [T6] . In Ref. [10] we have 
found that such an instability leads to a first-order phase transition even at finite temperature 
T between two phases of the dislocation - smooth and rough. Such transition is characterized by 
hysteretic behavior once dislocation length exceeds a certain threshold Lh growing with T |10j . 
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Figure 1. (Color online) Superfluid compress- 
ibility of the dislocation core. The resonant-type 
peaks occur at the thresholds F = F c (L,n), n = 
1, 2, 3 for n jog-antijog pairs creation. The parame- 
ters are L = 56, T = 0.05, up = 3.0. 
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Figure 2. Renormalized superfluid stiffness p s (T) 
and the velocity V s (T, F) of first sound normalized 
by their respective low-T values for different F 
(shown on the inset), L = 30, up = 3.0. Inset: 
the region of the dip (cf. Fig. 4 of Ref.[9]) showing 
its shifting with F. 
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Figure 3. (Color online) The left panel: power law dependence of the critical stress for several 
temperatures T = 0.006 — 0.013 (in the dimensionless units). The right panel: the same 
dependencies for all temperatures collapsed to a single master curve. The parameters Fq,Lq 
are shown for two different values of the Peierls potential strength. Error bars are shown for all 
points. Their bigger values on the right panel are due errors in values of 7. 



At sizes L < Lh the roughening is characterized the resonant-type increase of the dislocation 



compressibility n, Fig. [TJ As a result, speed of the first sound along the dislocation core exhibits 
a strong suppression shown on Fig. [2] for the first resonant peak. As can be seen, the strongest 
peak in Fig. Q] corresponds to the single pair creation. Thus, in what follows we will be ignoring 
the higher order processes. 



3.1. The critical stress for the dip 

The resonant-type increase of the compressibility occurs at some value of the critical stress 
F = F c which depends strongly on dislocation length and less strongly on temperature. With 
good accuracy such dependence has been found to be 

F C (L,T) = F {^j (6) 

where Fq, Lq are constants independent of temperature, and j(T) varies from unity in the 
limit of T = to 7 ~ 1.8 at temperatures approaching the temperature comparable to the 
energy of a jog-antijog pair A. Fig. [3] demonstrates such power law dependence for several 
temperatures. The constants Fq, Lq depend on the strength of the Peierls potential up, that is, 
on the parameters of the jog-antijog pair. In terms of the dimensional values, Fq is determined 
by a typical Peierls stress up ~ 10 — lOObar ( Fq ~ bap) and Lq is given by a typical jog-size. 
In what follows we will be considering stresses F which are much smaller than Fq. 



3.2. The dip temperature 

In the case when the free segment length is determined by the mean distance between the 
impurities as presented in Eq.(|2j), the condition of the resonance ([6]) can be eventually met for 
any (small) value of the external stress F. As will be shown below the temperature Td at which 
this occurs depends weakly (logarithmically) on the actual value of F. Indeed, substituting 
Eq.([2|) into Eq.© and solving for T in the approximation 7 ~ const we find 

E 

Td = "in [L (F /FyhX 3 -x d ]- (7) 

This equation is valid when < Lq{Fq/ F) 1 ^ X% — x^, < 1. Given « 1, this translates into 
the wide range of the stresses 

(A3L0) 7 < F/Fq < (X 3 L /x d y. (8) 

Since A3 /x^ ~ 1, Lq > 1 and F « Fq, it is reasonable to ignore x& in the denominator of 
Eq.([7j). Then, the dip temperature becomes much less than the activation energy 

Td = "ln [L (F /FyhX 3 ] ' (9) 

provided In Lq(Fq/ F) 1 / 1 A3 « — 1. It is important that T^ depends only logarithmically on 

3 He concentration and the values of the random stresses. As an example, taking A3 = 10 -7 and 
using F/Fq = W' 4 , L = 10, 7 « 1.5, E a = 0.5K, we find T d w 0.07K, with the log-factor being 
about 8. 3 He concentration could be varied in a wide range of values: from 10 -12 — 10 -15 in, 
practically, 3 He free samples to 10 -4 — 10~ 5 . Thus, despite the logarithmic dependence, a quite 
significant shift in the position of the dip can be observed. We find important studying such 
dependence experimentally. For example, for A3 ~ 10~ 12 and presuming the same protocol for 
preparing the sample, so that, the internal stresses remain essentially the same, one finds the 
dip temperature reduced by more than a factor of two. 



As discussed in Ref.[10J, the dip width is exponentially decreasing with the dislocation length 
(and eventually the resonant behavior turns into the hysteretic one). This conclusion is valid for 
a single dislocation. In the case of the dislocation forest with spatial fluctuations of the frozen 
in stresses the dip depth should rather be determined by a sort of inhomogeneous broadening 
- controlled by variations of F over a sample. If the average (mean square) fluctuation of the 
stress 5F is smaller than the mean square stress (F) itself, the estimate of the depth of the 
inhomogeneously broadened dip becomes 

5Td= 7 |ln [L (F /(F)yhX 3 ]\jF)- {W) 

Given that for practical values of the parameters the logarithm is of the order of 10, we conclude 
that 5T d «T d even if 5F ~ (F). 



4. Discussion 

Here we have proposed that 3 He impurities in combination with random frozen in stresses are 
controlling the dip-anomaly in the superflow rate through solid 4 He observed in Ref.|9j. For low 
density of dislocations the length of the dislocation segments L undergoing the stress induced 
roughening |10| is controlled by the inter-impurity distance along the dislocation cores rather 
than by the cross-pinning of dislocations. Due to the frozen in spatial fluctuations of stresses, 
there should generically always be segments meeting the condition for the resonant-type jog- 
antijog pair creation leading to a strong suppression of the flow. The exponential dependence 
of L on temperature insures that such resonances occur in a narrow temperature range (see 
Eq. (|10p ) around the logarithmically suppressed temperature Q. Thus, in contrast to a single 
dislocation in ideal crystal, in real samples the dependence on the average stress becomes much 
weaker. This may be the reason for a very weak dependence of the flux rate in the dip region 
found in Ref.|17j. Fig. 20. Under these circumstances we find important studying experimentally 
the dependence of the dip temperature on 3 He concentration. 
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